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Abstract 

We introduce the intuitive method to select an analytic Abel function of an ana- 
. lytic function f at a non-fixpoint. Due to the complexity of this method by involving 

^ ' matrix inversion of increasing size there is little known about its convergence. 

We show its convergence in the simplest but still complicated case f(x)=bx. 
We show that the obtained Abel function is, as expected, the logarithm to base b, 
independent on its development point. As a by-product we obtain a new polynomial 
approximation sequence for the logarithm to base b. 
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Q ■ 1 Introduction 

In the context of discussions of non-integer iterates of the exponential function there 
emerged a method which — in case of success — selects an analytic Abel function at a 
non-fixpoint. 

For a function / : G-^Cwe call a function a : G R f^ 1 (G) — > C Abel function of / iff 
^ . it satisfies the Abel equation 

a(f(z)) = a(z) + l (1) 



on its domain. Abel functions are an essential tool for non- integer/continuous iteration. 

>ri 

fM(z) = a-\t + a(z)), 



For a being bijective in an appropriate way one can define iterates by 



^ . they satisfy 

b ! 

. 5^ ! /M = / = /W o /W. 

for t being contained in some additive semigroup of C containing 1. Particularly p n ' is 
the n-times iteration/composition of the function / for positive integers n. 

We always consider Abel functions up to an additive constant, as one can see that 
if a — cti satisfies ([TJ then also ce(z) = cei(z) + c satisfies (pQ). However even up to an 
additive constant analytic Abel functions are not uniquely determined: If 9 is an analytic 
1-periodic function then a(z) = 9{z) + 9{a\{z)) also satisfies the Abel equation ([1]) which 
is easy to verify. 

There is an exhaustive theory about existence and uniqueness of analytic iterations 
(and the corresponding analytic Abel functions) developed at a fixpoint of /, see e.g. 
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Szekeres Ecalle jij or the monograph We refer to this method as regular iteration 
following Szekeres. 

As the exponential function e x has no real fixpoint, regular iteration is not applicable 
and quite different methods emerged aimed at obtaining real-analytic Abel functions 
anyway 0, EH, 0] • "Methods" here includes recipes with unverified outcome. 



For example some years ago Peter Walker [12] was proposing a way to calculate the 
powerseries of an Abel function a of the exponential f(x) = e x by solving an infinite linear 
equation system. His method was independently rediscovered in the lay-mathematical 
community (Andrew Robbins jsj), which documents a great interest for these kind of 
questions. 

His method works as follows: We consider the Abel equation 

a o f = 1 + a 

with formal powerseries a and / (in the hope that we obtain a with /(0) inside its 
convergence disk). We write the coefficient of x n in the formal powerseries / as f n . Then 
the formula for powerseries composition is (a o f) m = J2'^=o a nf n m where f n m is the 
coefficient of x m in the n-th power of /. The Abel equation can then be written as the 
infinite equation system in the coefficients of a: 

oo 
n=0 

If we subtract a m on each line m then we get the standard form of an infinite linear 
equation system: 

oo 

V (f n m - J m>n ) a n = J Wi0 m > (2) 




n=l 



The first column A m ^ is always as / = 1 and f° m = for m > 1. That's why we start 
the sum (and a n ) with n = 1 (and we know anyway that a may be determined merely 
up to CKq). 

Still this equation system must have infinitely many solutions, if any, of a containing 
/(0) in its convergence disk as we explained before. The intuitive method to solve this 
equation system — and hence to select one of the infinitely many solutions — is to solve 
the to N x N truncated equation systems 

N 

V (f n m ~ Im,n) OL™ = I mfi m = 0, . . . , N - 1 (3) 

n=l . 

-^T-m,n. 

for increasing N with the solution a[ N \ . . . in the hope that a n := limjv->oo a ^ 

exists for every n and the so obtained coefficient sequence (or formal powerseries) a 
is a solution of the untruncated equation system (or Abel equation) and has non-zero 
convergence radius. 
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Definition 1 (fjv, intuitive). For an infinite linear equation system Ax = b we denote 
the to N rows and N columns truncated matrix with A\n, similarly we denote truncated 
vectors. 

We call (x m ) m( zf$ the intuitive solution if the limit x m := limjv> m ,jv-»-oo (^Hv & Tjv ) 
exists for each m. 

We call (j2J) the (infinite linear) Abel equation system of / (developed at 0). We call its 
intuitive solution with «o = the intuitive (formal) Abel powerseries if existing. If a is 
analytic at (i.e. having non-zero convergence radius) we call the corresponding analytic 
function just the intuitive Abel function, written as .A^/]. 

We call the intuitive Abel function of g(x) := f(x+s)—s (for s inside the convergence 
disk of /) the intuitive Abel function of / developed at s, written as *4.J[/]. We call the 
function (3(x — s) the s-intuitive Abel function of /. 

The last part perhaps needs some explanation. Generally if is an Abel function of 
a conjugation g = h^ 1 o f oh then a = (3 o h~ x is an Abel function of /: 

ao f = fioh- 1 o f = /3 o goh- 1 = (1 + ^oh- 1 = l+a. (4) 

Particularly this is true for shift conjugations g(x) = f(x + s) — s. In this case the above 
(3 is the intuitive Abel function developed at s. By definition (3(s) = 0. The above a is 
the s-intuitive Abel function of /. 

Several until now unanswered questions arise here: For which / does the intuitive Abel 
powerseries exist (i.e. the coefficients converge)? Is A s independent on s in the sense that 
A x [f](x) — A x [f](x + s) is constant in xl Is it generally invariant under conjugation, i.e. 
for which h is A x [h~ 1 o f oh]— A x [f] o h constant? How does it relate to regular iteration 
at a (nearby) fixpoint? 

Besides the above questions there also arises the question whether this procedure gives 
the expected results for known elementary Abel functions of /. The most basic example 
being f(x) = bx with the Abel function a(x) = log 6 (x). 

2 Intuitive Abel function of f(x)=bx 

/ has already the fixpoint and it should be noted that the regular Abel function de- 
veloped at this fixpoint is indeed log b . The intuitive Abel function can however not be 
directly developed at fixpoint 0, because in this case the first line of our equation system 
is: Octi + 0a 2 + ■ ■ ■ — 1. 

So we proceed by calculating the intuitive Abel function (3 developed at s ^ 0, i.e. 
the intuitive Abel function of the shift conjugation g(x) = b ■ (x + s) — s which gives the 
s-intuitive Abel function L^ s (x) = f3(x — s). We will later see that L& )S is independent 
on s up to an additive constant. 

2.1 Solving the truncated linear equation system with a recur- 
rence 

In this subsection we solve the truncated equation system (J^D for / being the above given 
shift- conjugation g obtaining the recursive formula (j^J). We call the solutions (3^ instead 
ofaW. 
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In order to determine the occuring g n m (which, recall, is the m-th coefficient of the 
n-th power of g) we calculate: 

g(x) n = (6x + s^_Q) n = (^jd n ~ k b k x k 

hence the Matrix A in §2§ is given by subtracting the identity matrix from the matrix 
given by 

(I d d 2 d 3 \ 
n \,n- mh m eo . m b 2db 3d 2 b 

m ) d b ^ BU= b 2 3db 2 

\0 b 3 ) 

(which is also called the Bell matrix (or the transpose of the Carleman matrix) of g, 
see [H) and then removing the first column. We have to solve the equation system 
A \ N (P{ N \ pW) T = (1, 0, . . . , 0) =: u, where 



A- m ,n ( m ) d b Im,r. 



<m,0 



n=l...N 



d d 2 T\ fp[ N) \ 

\ :; 16-1 2db 3d 2 b 
b 2 -l 3db 2 



v/TV 




m=0...JV-l 



Then we equivalently change the equation system by multiplying each row n with s 7 



A' 



n 



in 



d n - m (sb) m -I m , n s r 



N=3 



d d 2 d 3 

sb - s 2d(sb) 3d 2 (sb) | pW 
(sb) 2 -s 2 3d(sb) 2 



As next step we equivalently change the equation system from A' to A"( N > (the N 
shall indicate that this matrix depends on the truncation size N while the previous 
steps where independent of N) by alternatingly adding up the lines onto the first line: 



A 



//(AO 



Yln=o(~ l) n ^rr The first line vanishes for 1 < m < N — 1: 



4 "(AO 



l)m s m + J2(~l) k ( ^ ) d m - k {sb) k 
k=0 ^ ' 

= ("!) 



= _(_i)™ s ™ + ( s (6_ l)-sb) 

d 

only the last entry at row m = N is non-zero: 

m— 1 



m+lgin _|_ 



A 



"(AO 
0,N 



E< 

n=0 



rn 



n 



d m - n {sb) k = {s{b - 1) - sb) rn - {-l) m {sbY 



s m (-l) m - 1 (6 m - 1) = (-l) w - 1 ((s6) 7V - I NtN s") 
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We multiply row with (— 1)^ 1 and rearrange the lines by moving row n + 1 one step 
up while row becomes the last row: 

'sb-s 2dsb 3d 2 sb \ 

v: 1 {sbf-s 2 3d(sb) 2 pW = 
(s6) 3 - s 3 J 

and drawing s m into (3, counting now rows and columns with first index 1: 




U„ 



n=l...N 



\n—mum 



b m - I 



h 



(AT) 



l m,N\ 



\N-1 



(b-l 2(6-1)6 3(6 - 1) 2 6\ ( sp[ N) \ 



N=3 









6 2 - 1 3(6 - 1)6 2 
6 3 - 1 



s¥ 2 N) 

3o(N) 




m=l...N 



The inverse matrix of a truncation of an upper triangular matrix is equal to the 
truncation of the inverse. We can give the n-th column of the inverse in terms of the 
values of the previous columns and the n-th column of the original matrix: 

1 ( ST 

\ k=m 

\ k=m / 

The multiplication with h^ N ' chooses the iV-th column of U~ l multiplied with sign 



\N-l. 



-1 



mo(n) 



) N 1 U m 1 N , yielding the recursion 

— (/^(-ir + E^'Q 1 

V k=m ^ ' 



1-6) 



n—k ik 



(5) 



2.2 A direct expression of the solution 

In this subsection we apply the technique of generating functions to obtain the direct 
(non- recursive) formula ([7j) for pm . Though I include the derivation, it is not necessary 
for the proof of the formula. So the uninterested reader may skip to proposition [1] where 
the actual verification of the formula takes place. 

We change from variable m to M as it will remain constant for our further considera- 
tions and the index m is needed to not run out of variables. Multiplying (jSJ) with 1 — 6™ 
and adding s M p^b n gives 

s M & = 4,m(-i) m + £ s M ^ (l) (1 - b)- k b k (6) 
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We further manipulate the equations 

>° nn _ t I i\Mun , un „M p{k)0 {± 0) 



v v ' k=M s v ' ^ ' 

to obtain the following recurrence in T n 



Tn _ T ( \\M , rp U 

7^-4,m(-1) +2^ t - (n _ fc)! 



fc=M 



of which we consider the generating function T(x) = J2^L T n x n . We get the left side of 
the last equation as the coefficients of T{x/b) 



OO ri~i 
J- n 



n=0 



and the right side is the multiplication of two formal powerseries (remember the formula 
(fg)n = ELo U-k9k), namely T and 



E 

leading us to 



;i - by 



X 3 = e (l-6)x 



3=0 J 



T(x/b) = (-x) M + T(x)e {1 - b)x . 
The reader may verify the following transformations for n > 0: 



i (1— Ma; 5^" 

T(b-^x) = (-xr IE" +T(x) e < 1 -^SS=o 6 



n _ i,-fc 



ytM 

,fe=0 



T(r^^) = (-*)" ( E * ^ ) + T(x)e^- (n+1) -V 



,k=0 

n e -xb- k ' 



-km-kM 



e-^~ (n+l) T{b-^x) - T(x)e- ta = (-a;) M £ „ w 

\fc=0 

= { _ x) m y {_^)y h - 

m=0 fc=0 
00 (_ X ^+ M £-(m+M)(n+l) _ j 



m=0 



e- fe:ES/ T(^) - T(x)e" ba; = 



°° ^_-Qm+M ytn+M _ X m+M 

6-(m+M) _ ! 
m=0 



E m , ( U + M)-l) 
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Letting now y = b n — > by n — > oo for \b\ > 1, considering S and T being continuous 
at and S(0) = T(0) = we get: 



oc 



T(x) = e bx S{x) = e bx 



(m-M)\(b- m - 1) 

m=M v ' v ' 

with the coefficients given by formal powerseries multiplication: 



fe=M 



(n-A;)!(A;-M)!(6- fc -l) 



n! 

k=M 



k \M 1 - b k 



So the direct formula is 



*?-<--i:(;)G)££ « 

fc=m V / \ / 

Proposition 1. T/ie direct expression ([7]) satisfies the recurrence ([5]) /or any 6 G C iaai 
zs not a root of unity (i.e. b n ^ 1 /or any n > 1), and is hence the solution of the truncated 
Abel equation system A\ n [3^ = u\ n of g(x) = b ■ (x + s) — s for s / 0. 

Proof. We just fill the direct expression ((7j) in equation fl6]) (which is equivalent to the 
recurrence (jSJ)) and show by equivalent transformation that the recurrence is satisfied: 

= LM-ir + e t $ ( k M ) ^ © (i - 

k=Mk'=M \ / \ / \ / 

1 )" + E©^E(l)(t)(-^ 

ui—li/r \ / fc_l./ \ / \ / 



'n,M< 



fc'=M v ' fc=fc' 
k'=M \ / \ / fc=fc' \ 7 



<n,Mv 



i ' m+ i:(S(m){ 

fc'=M \ / \ / 



n\fk'\ {-bf 



h k> 



Then we replace s M (3^ on the left and subtracting the right sum 

fn\fk'\ (-rf -(-bf 



k'=M 



k')\MJ l-b k ' 

fc'=M \ / \ / 

;)e(;::)h^u (-D 

_l)M Q iV- m = 7 (-1 



M 



M 



M 
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And this is indeed a true statement, considering fc = for k > 1 and 0° = 1. □ 

2.3 Convergence of the polynomial approximation 

Now, that we obtained the truncated solutions /3^ N \ we want to see whether the limit 
(3 n = liHiAr^oo /% (according to our definition [1] of "intuitive solution") exists for each 
n, which would be then the n-th coefficient of the intuitive Abel function (5 of g; where 
a(x) = (3(x — s) would be the s-intuitive Abel function of f(x) = bx (compare fl4])), which 
we want to prove to be a(x) = log fe (x) + c for some c possibly depending on s. 

The coefficient wise convergence would be a direct consequence of the pointwise con- 
vergence of the polynomial approximations a^ n \x) = Ylm=i@™\ x ~ s ) m ^° a ( x )- I n 
the following subsection we prove this convergence by showing that 
converges to log b (a;). Our efforts culminate in the summarizing theorem [TJ 



a {n) (x) 



£(— >-™£(:)Ct)S 

m=l k=m x ' v ' 

e(;)£$i>/.-0 

k=l v 7 m=l v 7 

' n \r i\k+i l ~ ( x / s ) k 



k=l 

which is the s-free function applied to x/s: 



b k 



a^(x) :=a( S x) = X;( r ")(-l) 
k=i ^ 7 



n\, ^ k+1 l-x k 



l-b k ' 



With little effort, 



E(?)(-irv=i-(i-y) B , (9) 
k=i ^ 7 

we can compute the value of a^(b m ) 

n r \ m—1 m—1 

aW(6 m ) = ^ ( J E ^ = E 1 " ^ ~ ^" 

fc=l ^ 7 t=0 i=0 

lim 5 (n) (6 m ) = m 

n—toc 

confirming our hypothesis that lim^oo = log b (x). However to prove it (for a 

series with f(a n ) = n for all n > 1 that is not the logarithm see e.g. Euler jij) we need 
to show that lim^oo aS n '{b x ) = x also for non-integer x. To be careful we restrict x and 
b from here throughout this section to < b < 1 and < x < 1. We have a look at the 
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series expansion of for < y < 1: 



1-2/ 



J + l 



(y - i) J 



E 

3=0 

oo 

7=0 VJ 7 i=0 



;-i) j -y 



(10) 

(ii) 

(12) 



Now substituting y = b k 



ky 



k=l \ ' 

And knowing that Y2=i (T)(-l) k+1 b ki = 1 - (1 - b l ) n we write 



(13) 



a 



(n) 



We split the sums into two parts at the first minus of 1 — (1 — b l ) n . Considering 
Ei=o 0(-l) j_i l = ° j > where ° j = for j > 1 and 0° = 1, we get 



«<»>(&") = * - E + x J E ( ■ J (-^(i - & T • 



R' 



(n) 



Obviously lim^oo i?^ = for j > 1 because |1 — 6 l | < 1 for each i. In the remainder of 
this section we show that the sequence (of sequences) converges not only point-wise 
but uniformly to (0, 0, ... ) in the supremum norm \ \v \ \ = sup JgN \vj\ which then implies 
that we can swap taking the limit in n with the limit in j. 



So we show that for each e > there is an uq such that R 
For j = 1 we find an n = n such that 

«r|2£(0ii-&r< £ 



(n) 



< e for all j and n > uq. 



Noticing that dj jTl is decreasing in the second index n the above equation is also valid for 
any n > n . 

By binomially expanding the power with exponent n and unexpand it into a power 
with exponent j we reformulate the expression of R^ to: 



R (n) = ^y-nS^ (f\ (-l)«- fe (l - b k Y 
k=0 ^ ' 
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and so obtain that 



R 



(n) 



which is now j. Hence 



R 



(n) 



Rn^ < d n j where d n j is decreasing in the second argument 
< e for all n > n and j > 1. 



Now, to finish the proof we show that 



lim ( 

n->oo ^ — ' V 7 + 1 



# = 



As prerequisite we need the well known 

Preliminary 1. Yl™=o O- 2 ™ conver 9 es absolutely (to (z + 1) K ) for all z with \z\ = 1 if 
ft(/c) > 0. 

Proof. Done via comparison with the Riemann zeta function by the inequality | (*) | < 
Corollary 1. JTie Kmi* t K := J2?=i Q+i) 



exists /or every k > 0. 



Proof. Follows from letting z — 1 in the preliminary. 



□ 



Then for a given e > choose no such that it- < e/t x for all j > 1 and n > no and 
obtain: 



(n) 



x 



— = £ for n > n Q 



So we have established that \im. n ^. 0O 6t( n \\y c ) = x for x > 0. Which has the consequence 
that a(y) := lim^oo <5f( n )(y) exists for all < y < 1 and is the inverse function of x (->■ fr 1 . 
We finally summarize all our findings: 

Theorem 1. TTie (unique) polynomial ft n \z) of degree n that satisfies (3^(0) = and 
the Abel equation 

(3 {n) (b ■(z + s)-s) = l + ft n \z) (s ± 0, b k ± Vik > 1) 

is given by (5^ n \z) = a^ n \z + s) = oS n \z/s + 1) where 

— (^V_i) fe+1 — 

\k p ' l-b k 

k=l v 7 

and for all x, b G (0, 1) we have the convergence: 



lim aS n \x) = log b (x). 



(14) 



a(x) = a(x/s) = log 6 (x) — log fe (s) is the s-intuitive Abel function of f(x) = bx. (It is 
independent on s up to an additive constant.) 
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3 Comments 



The most urgent questions to develop the mathematics of the intuitive method are already 
listed in the introduction. Here only some side notes: 

Numerically it appears that convergence of at is also achieved for \x/b — 1| < 1 in the 
case b > 1 which points towards possible improvements of the theorem. 

The more interesting question about the convergence of the approximating polynomi- 
als of the intuitive Abel function of f(x) = e x seems out of reach to solve with these rather 
elementary techniques. Numerically at least it seems that the coefficients do not converge 
uniformly but have a point-wise limit which is invariant under shift conjugations. 

The author thanks for the stimulating discussions on [10]. Without them this paper 
would not have come into existence. 
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